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The electronic structure and magnetism in the sodium nickelate NaNi02 in the low-temperature 
phase is studied from density-functional calculations using the linear muffin-tin orbitals method. 
An antiferromagnetic solution with a magnetic moment of 0.7 /is per Ni ion is found. A ferrodis- 
torsive orbital ordering is shown to occur due to the Jahn- Teller distortion around the Ni +3 ion 
in agreement with the orbital ordering inferred from neutron diffraction. While the intralayer ex- 
change is ferromagnetic, the interlayer exchange is weakly antiferromagnetic, mediated by a long 
Ni-O-Na-O-Ni superexchange path. 



Since their discovery in the early fifties^, the nickelates 
have been the subject of attention due to their many ap- 
plications such as the base materials in batteries. More 
recently, the family of compounds L^Mei-aNiC^ (Me 
being a metal) have attracted special interest due to their 
triangular lattice structure which makes them good can- 
didates as frustrated magnetic systems. In fact, magnetic 
or orbital frustration mechanisms are usually invoked in 
explaining the unusual behavior of LiNiC>2- Many differ- 
ent scenarios have been proposed to explain the absence 
of orbital and magnetic orderings in LiNiC>2, such as spin- 
glass^ quantum disordered statej&a spin-orbital liquid, 5 
impurity effect, 6 or frustrated anti-ferromagneti In com- 
parison, NaNiC>2 seems to show none of the unusual be- 
haviors of its sister compound including the magnetic 
behavior For example, while in LiNi02, no long-range 
magnetic order has been observed, NaNiC>2 in contrast is 
a type A antiferromagnet (ferromagnetic layers stacked 
on top of one another and coupled antiferromagnetically) . 
It has been a long puzzle as to why their magnetic prop- 
erties are so different, in spite of the fact that the two 
compounds are very similar. 

In this paper, we study the electronic structure of 
NaNi02 using density-functional methods, with the goal 
of gaining insight into the physics of this family of mate- 
rials. The nature of orbital ordering as well as the origin 
of the magnetic exchange are also discussed. 

The high-temperature hexagonal crystal structure of 
NaNi02 (space group R3m, no. 166) is shown in Fig. 
^ The compound undergoes a structural transition to 
a lower-symmetry monoclinic structure with the para- 
magnetic space group C2/m (no. 12) at about 450 K, 
below which the oxygen octahedra become elongated. 
The magnetic transition is at a much lower temperature 
Tyv = 20K 9 , below which the A- type antiferromagnetic 
structure occurs. The layered structure may be viewed 
as an arrangement of slightly elongated NiC>6 octahedra 
separated by Na sheets. In the low-temperature struc- 
ture, the Jahn- Teller (JT) distortion leads to different 
Ni-0 bond lengths: four short bonds of 1.91 Aand two 
long ones of 2.14 A. The lattice parameters are listed in 
Table [|] In the low-temperature phase, NaNiC>2 shows 
the ferro-distorsive orbital ordering caused by the Jahn- 
Teller distortions of the NiC>6 octahedra, which disap- 




FIG. 1: High temperature crystal structure of NaNi02- Cal- 
culations reported here are for the low-temperature phase 
(C2/m), but the high-temperature structure (R3m) is shown 
here because it is simpler to visualize. The former is obtained 
by distorting the latter, including a Jahn- Teller stretching 
of the NiC>6 octahedra. The magnetic ordering is antiferro- 
magnetic type A, with ferromagnetic interaction within the 
layer and antiferromagnetic interaction between the layers. 
The line connecting the Ni-O-Na-O-Ni atoms indicates the su- 
perexchange path and the rectangle joining the four Ni atoms 
has reference to the plane of the contour plot in Fig. 4. 



pears in the high-temperature phase. 

The band-structure calculations were performed for 
the low temperature structure using the local spin- 
density approximation (LSDA) to the density functional 
theory (DFT) as well as the "LDA+U" approximation, 
using the linear muffin-tin orbitals (LMTO) methodi2*ii. 
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FIG. 2: Density-functional energy bands for the antiferro- 
magnetic NaNi02 in the low-temperature structure obtained 
from the "LDA+U" calculations. The unit cell has two for- 
mula units: ( NaNiO^^- The e 9 bands are split near the 
Fermi level due to the Jahn- Teller and exchange interactions. 



For the antiferromagnetic structure, the symmetry is not 
reduced further and the magnetic unit cell is also mon- 
oclinic with space group C2/m with two formula units 
per unit cell. Within the LMTO atomic sphere approx- 
imation (LMTO-ASA), the antiferromagnetic (AF) type 
A structure was found to be the ground state, lower in 
energy than both the ferromagnetic (FM) and paramag- 
netic structures. 

In the "LDA+U" method^ the Coulomb interac- 
tions between the localized electrons (Ni(d)) are sub- 
tracted from the LDA total energy and a "more proper" 
Hubbard-like term is added, so that the corrected energy 
functional reads: 



E = E 



LDA 



UN(N- 1/2) 



i¥=3 



(1) 



where N is the total number of the localized electrons, 
rii is the occupancy of the i th localized orbital, and U is 
the screened Coulomb energy. The orbital energies are 
obtained by taking the derivatives with respect to the 
orbital occupation, so that 



dE/dn i = e LDA + U{l/2-n i ). 



(2) 



Thus with respect to the LDA results, the energies of the 
occupied orbitals are shifted by —U/2, while those of the 
unoccupied states are shifted up by the amount U/2. In 
the actual calculations, one works with a slightly more 
complicated functional that includes properly the direct 
and the exchange Coulomb interactions, parametrized by 
the screened interactions U and J. 

The band structure of NaNiC>2 in the AF-type A mag- 
netic structure calculated with the "LDA+U" method, 



TABLE I: Atomic positions and sphere radii used in the 
LMTO calculations. Lattice parameters for the monoclinic 
unit cell (space group C2/m) are a — 5.311 A, b = 2.840 A, 
c = 5.568A, and j3 = 110.44° (Ref£). 
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with U — 5 eV and J = 0.5 eV, are shown in Fig. 
While the LSDA results showed a metallic band struc- 
ture because of a slight overlap between the bands de- 
rived from the e g \ and e g J. Ni(d) states at the same Ni 
site, the overlap disappears when the correlation effects 
are included within the "LDA+U" formalism. Only a 
small value of U is needed to open up the gap, suggest- 
ing a robust insulating gap in the compound resulting 
from electronic correlation. 

The bands indicate a low-spin state , with the tig states 
being completely occupied while the e g states are only 
quarter-filled (i|_ej). As indicated from the densities- 
of-states (Fig. the t 2g and e g bands are split by a 
strong crystal-field splitting, while the ej and e g are split 
by the exchange coupling by the amount A ex « 0.7 eV. 
The Ni(eg) occupancy makes the system JT-active in the 
low-temperature phase, lifting the degeneracy of the e g 
levels, with a JT energy of Ajt ~0.6 eV. 
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FIG. 3: Electronic Density of states for NaNi02. Ajt and 
A ex refer to the Jahn- Teller and the exchange splittings, re- 
spectively. 

The calculated magnetic moment is fi =0.51 /is/Mn 
obtained from the LSDA, which increases to 0.71 /x^/Mn 
in the "LDA+U" results, as in the latter, the slight over- 
lap between the e g f and e g J, disappears, increasing the 
occupancy of the spin f states. This is lower than the 
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FIG. 4: (Color online) Charge density contours of the occu- 
pied Ni (d) bands (about 0.8 eV wide occurring just below 
Ef) showing the ferro-orbital ordering of the Ni(e 9 ) orbitals, 
as obtained from the LSDA calculations. The plane of the 
plot is the diagonal plane, defined by the rectangle drawn 
here and in Fig. Arrows indicate the spins of the Ni 3+ ions 
and the line joining the Ni-O-Na-O-Ni indicates the superex- 
change path. Contour values are given by p„ = po X 10 -ni5 , 
with po = 0.708 e~/ag, 5 = 0.575, and n = 0, 7. 



nominal magnetic moment of 1 ps/Ni because of the 
strong hybridization between the 0(2p) and the Ni(3d) 
states. 

We have studied the orbital ordering by examining the 
occupancy of the various Mn (d) orbitals. In a frame 
of reference where the z axis points along the long Ni-0 
bond, we found the occupied Ni(d) electrons to be pre- 



dominantly of 3z 



character. The charge density 



contours for the occupied Mn (d) states, plotted in Fig. 

shows the ferro-orbital ordering of the 3z 2 — r 2 or- 
bitals, with significant hybridization coming from the O 
(2p) states. Fig. 5 shows a sketch of the orbital order 
as obtained from the density-functional charge density. 
The charge density for the undistorted high temperature 
structure was also computed and we found no orbital or- 
der for that case. 

The magnetic exchange interactions, the ferromagnetic 
Jp within the layer and the antiferromagnetic Jaf be- 
tween the layers are both weak, with the measured values 
being Jp = 13 K and Jaf = —1 K*^ The exchange inter- 
actions are weak because the mtralayer exchange involves 
a 90° Ni-O-Ni bond, while for the mierlayer coupling, it is 
the long Ni-O-Na-O-Ni superexchange path that makes 
the exchange weak. The 90° Ni-O-Ni exchange within 
the layer has been examined in great detail in two recent 

For the exchange between the layers, we treat the Ni- 
O-Na hopping as an effective Ni-Na hopping between the 



Ni(3z 2 



2 ) and Na(s) orbitals on the Ni-O-Na-O-Ni 




FIG. 5: Sketch of the orbital order in NaNiCh extracted from 
Fig. 4. The shaded lobes indicate a positive wave function, 
indicating the anti-bonding interaction between the Ni (3z 2 — 
r 2 ) and the 0(p z ) orbitals. 



fective Ni-Na-Ni path then reads 



E 



Sin-, 



Uin^n^, (3) 



where t is the effective hopping between the Ni and 
the Na atoms, £, is the on-site energy, and the on-site 
Coulomb energies are Ud and U s for the Ni (d) and the 
Na (s) orbitals, respectively. 

The exchange interaction Jaf between Ni atoms is ob- 
tained by calculating the ground state of the Hamiltonian 
Eq. |3 for the ferromagnetic (FM) and antiferromagnetic 
(AF) configurations of the Ni (e^) spins and then taking 
the difference: Jaf — -E^l — -^TT- For the FM case, we 
have two spin up electrons occupying the Na (s) and the 

Ni (e^) states leading to a three-dimensional configura- 
tion space. Similarly, for the AF case, we have a 9 x 9 
Hamiltonian (one spin up and one spin down electrons 
distributed among three orbitals of either spin, so that 
3 Ci x 3 Ci = 9). The ground-state energies may be ob- 
tained by diagonalization of the Hamiltonians or by using 
the fourth order perturbation theory. The latter yields 



Jai 



2t 4 



U s + 2A 



sd 



i 

~Ua 



(4) 



superexchange path. The model Hamiltonian for the ef- 



where A S( j is the charge transfer energy from Ni (d) to 
Na (s). The first term is due to the exchange via the in- 
termediate state with the two electrons on the Na atom, 
while the second is due to hopping between the two Ni or- 
bitals via the Na atom. The exchange coupling is clearly 
antiferromagnetic. Taking typical values: t w 0.3 eV, 
A s d » 5 eV, Ud ~ 5 eV, and U s « 2 eV, the magnitude 
of the exchange is Jaf ~ —0.1 meV, which is the right or- 
der of magnitude as compared to the experimental value 
of -1 K. 
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